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New solutions are obtained for the planar motion of three equal mass particles on a common periodic 
orbit with finite total angular momentum, under the action of attractive pairwise forces of the form 
l/r p+1 . It is shown that for —2 < p < 0, Lagrange's 1772 circular solution is the limiting case of 
a complex symmetric orbit. The evolution of this orbit and another recently discovered one in the 
shape of a figure eight is investigated for a range of angular momenta by a novel numerical method. 
Extensions to n equal mass particles and to three particles of different masses are also discussed 
briefly. 
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Introduction 

Recently C. Moore! applied the theory of braids to find 
a classification of periodic orbits for the motion of several 
particles moving in a plane under the action of pairwise 
attractive power-law forces, obtaining numerical results 
for the case of 2 and 3 strands. In particular, for three 
particles of equal mass, he found a solution where these 
particles travel on a common orbit in the shape of a fig- 
ure eight. This orbit was found-also independently by A. 
Chenciner and R. MontgomeryB,u who have proved its ex- 
istence for the case of zero total angular momentum and 
inverse square forces. Moore's numerical result that this 
orbit is stable has also been confirmed by computations 
of C Simoo. Previously, one of the simplest known peri- 
odic orbits without collisions for this problemtj was La- 
grange's classic solutioro where the three particles form 
the vertices of an equilateral triangle with each particle 
moving on an elliptic orbit with focus in the common 
center of mass. In this paper we show that for attractive 
power-law forces of the form l/r p+1 with — 2 < p < 0, 
the circular Lagrange orbit is the limit of a family of 
orbits which are periodic in a uniformly rotating frame 
of reference. For small deviations, we find analytically 
that this orbit is a distorted ellipse which, as the angular 
momentum decreases, resembles a distorted figure eight 
shape, but it does not approach the recently discovered 
zero-angular-momentum figure eight orbit which has a 
different winding number. For this figure eight orbit, 
which exists for — 2 < p, we obtain an analytic approxi- 
mation, and we prove by perturbation theory that it can 
be extended to finite angular momentum. We also inves- 
tigate the evolution of these two orbits for a finite range 
of total angular momentum by a novel iterative method, 
and consider briefly the case of such periodic orbits for n 
equal mass particles, and for three particles of different 
masses. 



Computational method for periodic orbits 

Our method is based on expanding in a Fourier series 
the Cartesian coordinates for the periodic orbit on which 
all the three equal mass particles travel, and then finding 
the extremun of the action for this system with respect 
to the Fourier coefficients in a frame rotating with a fixed 
frequency tt. We obtain analytic approximations by per- 
turbation theory, and numerical solutions by a novel it- 
erative method which converges rapidly. Setting 9 = cut, 
where t is the time, oj — 2t:/t is the frequency, and r is 
the period of the orbit), we have 

k—n 

x(9) = a s (k )sin(k 9) + a c (k e )cos(k e 9) (1) 
fe=i 

and 

k—n 

y{9) = b s (k e )sin(k e 9) + b c (k )cos(k 9), (2) 

where k Q = 2k — 1 and k e = 2k, excluding multiples 
of the integer 3. In practice, the number of terms n in 
the Fourier sums, Eqs. [l] and |], depends on the desired 
accuracy of the solution. The Cartesian coordinates of 
each of the three particles is then given by Xi = x(9i) 
and yi — y(9i), for i = 1,2,3, and the invariance of the 
action under the cyclic transformation 1 — > 2 — > 3 — > 1 
is satisfied by the requirement 9\ — 9, 9^ = 9 + 2ir/3 
and #3 = 9 + 4ir/3. These conditions are satisfied by the 
equations which determine the extremun of the action in 
a frame rotating with frequency Q in the plane of the 
orbitCI which takes the form 

,2tt 

A = / d9[K(9)-P(9) + -n 2 I(9) + nL(9)}, (3) 
Jo 2 

where K(9) is the kinetic energy, P(9) the potential en- 
ergy, 1(9) the moment of inertia, and L(9) the angular 
momentum in the rotating frame. These are given by the 
following expressions: 
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^) = -A(i + 4 + 4), (5) 
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where r.y = ^/ (x* — Xj) 2 + (j/i — j/j) 2 is the distance be- 
tween particles i and j, and p is the power of the force 
law, 

The extrema of the action, Eq. [|, are obtained by find- 
ing the Fourier coefficients for which the partial deriva- 
tives of the action with respect to these coefficients van- 
ish, 



dA 



da s (k a ) 



dA 



db c (k ) 



u(k )a s {k o ) - f s (k a ) - v(k )b c (k o ) = 0, (8) 



u(k )b c (k o ) - g c {k ) ~ v(k )a s (k o ) = 0, (9) 



dA 



db s {k e ) 



= u(k )b s {k e ) - g s {k e ) + v(k )a c (k e ) = 0, (10) 



We obtain numerical solutions of these equations by 
starting with the Fourier coefficients to some conjectured 
approximate n period orbit, e.g. according to Moore's 
classificationEI, and then evaluating the partial deriva- 
tives of the action, Eqs. Bjlll If the extrema of the 
action is at a maximum (minimum), an improved orbit 
is then obtained by changing each coefficient in propor- 
tion to plus (minus) the corresponding partial derivative 
of the action. Thus, starting with some initial value for 
the coefficient a s (k ), we obtain a new value a s {k ) by 
the relation 



a' s (ko) = a s (k ) + 



Ss dA 
u(k ) da s (k ) 



(16) 



with similar expressions for the other Fourier coefficients 
where Ss is a positive (negative) parameter if the action is 
a maximum (minimum) . This procedure is then iterated 
until the partial derivatives of the action are sufficiently 
small to any desired accuracy. Our method differs in 
an essential way from the standard method of steepest 
descent by the important factor l/u{k) which is essential 
for the convergence of the iterations with moderate values 
of 5s (order .1). 

To obtain analytic approximations of Eqs. p[]Tl|, it is 
convenient to determine the frequencies lo and Q by fixing 
the values of a s (l) and b c (l), which give the scale of the 
orbit and the magnitude of the total angular momentum 
in the rest frame, Lt = L + QI, where L is the mean 
angular momentum in the rotating frame, Eq. |t], 
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fc=l 



-fc a s (fc )6 c (fc ) + k e a c (k e )b s (k e 



and I is the mean moment of inertia, Eq. O, 



(17) 



dA 



da c (k e ) 



u(k )a c {k e ) - f c {k e ) + v(k )b s (k e ) = 0, (11) 



where u(k) = k 2 uj 2 + Q 2 and v(k) = 2kujQ. Here the 
functions 



fs{ko) = 



Sc(k e 
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d6— — cos(k e 9), 
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dP 



g s (k e ) = — d0— — sin{k e 6) 



9c(k ) = 
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depend also on all the Fourier coefficients. 
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I = | E a s (k ) 2 + a c (k e ) 2 + b s (k e f + b c (k ) 2 . (18) 
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We have 
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where 



and 



/*(!)- g c (l) 
(l-&c(l)) 

/«(!) + g e (l) 

(i + Mi)) 



(19) 
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are obtained from Eqs. |8| and ^ for k Q — 1 and a s (l) — 
1. These frequencies scale with the size s of the orbit 
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according to the generalized Kepler relations that oj 2 s p+2 
and Sl 2 s p+2 are constants. n 

The generalization of Lagrange's solutions is obtained 
by setting the Fourier coefficients a c (k e ) and b s (k e ) equal 
to zero, which according to Eqs. [ll] and implies that 
f c (k e ) and g s (k e ) also vanish. Small deviations from 
Lagrange's circular solution can be obtained by setting 
a s (l) = 1 and b c (l) = 1—6, and expanding Eqs. |^ and 
in a power series 5. In the limit of vanishing 6 



3P/4 



(23) 



which shows that these solutions exist only for p < 0. In 
particular, for the special case that p — —2 corresponding 
to harmonic forces, ui = V3, f2 = 0, and all the Fourier 
coefficients a s (ko) and 6 c (fco) vanish for fco > 1. 



Numerical solutions and analytic approximations 

Two numerical solutions of the extended Lagrange 
problem are shown in Figs. 1 and 2, where segments 
of the orbit occupied by each particle during a third of 
a period are indicated by a dark line, a dotted line and 
a dashed line respectively, and ten successive locations 
on these segments are marked at equal intervals of time 
to illustrate the relative position and variable speed of 
the three particles moving on this orbit. In this solution 
the angular momentum of the orbit decreases as b c (l) 
decreases, and the approximate elliptical orbit starts to 
get pinched in the middle. This is shown in Fig.l for 
p = — .05 at the special value 6 C (1) = .049 where the 
orbit also develops cusps at the two ends. At these val- 
ues of the parameters each particle comes momentarily 
to rest at these points because the attractive forces due 
to the other two particles balance the centrifugal force in 
the rotating frame. 

Further decrease of the value of b c (l), corresponding 
to decreasing angular momentum, requires that the y- 
component of velocity at these ends change direction, 
which is possible only if an additional lobe appears at 
each end. A further decrease of 6 C (1) leads to an addi- 
tional lobe at the center, as shown in Fig. 2 for p = —.05 
and 6 C (1) = .0171. 

The recently discovered] -u zero angular momentum 
periodic orbit in the shape of a figure eight is a solution 
of Eqs. H and |l0| with the Fourier coefficients a c (k e ) and 
b c {k ) set equal to zero. According to Eqs. lTj and ||, 
this implies that f c {k e ) — g c (k ) = 0. For the special 
case that p = 2, it can be shown that the moment of 
inertia, Eq. |^, is a constant. Consequently the shape 
of this orbit can be determined from purely geometrical 
considerations by imposing the condition that the Fourier 
coefficients of an expansion of this moment and that of 
the angular momentum vanish. Keeping terms in the 
sums in Eqs. [j] and [2] up to n=4, we obtain 




FIG. 1. An extension of Lagrange's circular orbit for three 
equal mass particles labeled here 1,2 and 3 , showing their 
location at ten equal time intervals during one third of a pe- 
riod. For example, when particle 3 reaches an extremun along 
the x-coordinate, its y-component of acceleration vanishes by 
symmetry, and therefore the location of the other two par- 
ticles, 1 and 2, lie on the vertices of an isosceles triangle as 
shown here for the case p — — .05 and fe c (l) = .049. 




FIG. 2. Extended Lagrange orbit with crossings near the 
center and additional lobes at the two ends, for p = —.05 and 
6 C (1) = .0171 
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6 s (4) = ^-|^(46 + 976 2 + 49& 4 ), 



o.(5) = ~ - ^j(25 + 496 2 + 2 26 4 ), 



a. (7) 



2d 2 ' 



6 3 

0.(8) - -g,, 



76 2 



(24) 



(25) 



(26) 



(27) 



where 6 — 6 S (2) and d = 5 
determined by setting cj = y/ f s ( 1 ) , Eq 
Eq. pi for jfe, 



The magnitude of & is 
|{| and solving 
2 which takes the form b = h(b), where 



h(b) 



9s(2) 
4/ s (l)' 



(28) 



We find that h(b) is a positive monotonically decreasing 
function for b > which diverges as b approaches zero, 
and vanishes as b increases indefinitely. Therefore, there 
exists only a single solution which we find at b = .3793 
corresponding to w = 1.448. This result is in good agree- 



and 
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which 



ment with our numerical solution of Eqs. 
gives b — .3776 and u> — 1.452. The resulting figure eight 
orbit obtained from the Fourier coefficients given by Eqs. 
p4| - ^ is shown in Fig. 3, which is graphically indis- 
tinguishable from the orbit obtained from the numerical 
solution of Eqs. || and [l0|. Surprisingly, these analytic 
coefficients also give a good approximation to the zero- 
angular-momentum figure eight orbit for p < 2, because 
it turns out that the,-moment of inertia is approximately 
constant in this caseo. As p decreases we find that b de- 
creases and vanishes when p — — 2, because in this limit 
of a harmonic force law the figure eight orbit becomes a 
linear collision orbit. 

The zero-angular-momentum orbit provides the basis 
for a perturbation calculation of the corresponding orbit 
for small values of the angular momentum. Expanding 
the coefficients a c {k e ) and b c {k Q ) in powers of b c (l) = b, 
these coefficients together with the unperturbed coeffi- 
cients a s (fc Q ) and b s (k e ) determine the angular momen- 
tum of the orbit. For example, for p — 1 we find that 
il = 1.6476, and substituting this value of and the 
lowest order Fourier coefficients in Eqs. [l7| and |lj we 
obtain L t — —.52766. Solutions for larger values of b are 
obtained numerically by iterating these equations, fixing 
a s (l) and starting with some arbitrary value for 6 S (2), 
while all other Fourier coefficients are initially set equal 
to zero. As b increases the figure eight orbit becomes 
more asymmetric with respect to the y-axis, as expected 
from the increasing action of the centrifugal and Coriolis 
forces in the rotating frame of reference, and at a critical 
value of the angular momentum the y-component of the 
velocity of the particle at the end of the orbit in a di- 
rection opposite to the rotation vanishes giving rise to a 




FIG. 3. Two analytic approximations to the figure eight 
periodic orbit with zero angular momentum for p = 2. The 
dashed line is the solution for n=l, a lemniscate, and the solid 
line is the solution given by Eqs. 23-26 for n=4 



cusp. A further increase of the angular momentum then 
requires that this component of velocity change sign at 
this point, which give rise to an additional loop in the 
orbit, as shown in Fig. 4. 

Finite angular momentum solutions can also be ob- 
tained by rotating the figure eight orbit about an axis in 
the plane of the original orbit leading to a three dimen- 
sional orbit as shown in Fig. 5 

Solutions of our equations provide initial conditions for 
a numerical integration of the equations of motion which 
can then be applied to study the stability of these peri- 
odic orbits with respect to small perturbations. A prelim- 
inary investigation indicates that some of these solutions 
are stable over .a restricted range of angular momentum 
and values of pS. Our method can readily be extended 
to investigate similar periodic motions for an arbitrary 
number of equal mass particles moving on the same or- 
bit, many of which have been discovered recentlya . For 
example, any odd number n a of particles can travel on 
a figure eight orbit similar to that shown in Fig. 3. We 
find that this orbit, Fig. 6, approaches a unique limit 
as n — > co, but the frequency increases logarithmically 
with n . 

In contrast, for the extension of Lagrange's circular 
orbit each additional particle leads to a different orbital 
curve. Finally we have also applied our method to study 
the case of three particles with different masses which 
must move, however, along different orbital curves. An 
example is shown in Fig. 7 which is an extension of an 
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FIG. 4. Asymmetric figure eight periodic orbit with finite 
angular momentum in a rotating frame for p = 1 and b — .38. 
The additional loop shown at the left end of the orbit forms 
when b > .34. 
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FIG. 5. Three dimensional figure eight orbit with finite 
angular momentum along the z-axis. 
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FIG. 6. Figure eight orbit with 201 particles 

orbit discovered by MooreEl for the case of equal mass 
particles. 
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FIG. 7. Periodic orbits for three particles with masses 
1.0, .5 and 1.5 at initial positions labeled 1, 2, and 3 respec- 
tively, under the action of inverse square pairwise forces. Par- 
ticles 1 and 2 move initially clockwise, while particle 3 moves 
anti-clockwise. Positions are shown at 1/16 of a period, and 
the vertex of the triangle shows their relative position after 
1/8 of a period. 
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